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Abstract 

Mass spectrum of O"*""*" glueballs is produced using a dual supergravity theory we proposed 
for pure M = 1 SU (N) gauge theory in four dimensions in the large limit in the IR. The 
glueball states are expressed in terms of Whittaker functions. The spectrum is discrete and a 
function whose roots give the masses is written. The ratios of the masses are independent of 
the parameters of the theory and comparison to recent non-supersymmetric large lattice 
QCD data available for the lowest three states shows agreement to within five percent. 
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1 Introduction 

One important development in theoretical physics in recent years has been the finding of 
examples of dual gauge and gravity theories following [U 121 E] in which a gauge theory 
in strongly coupled nonperturbative region has a weakly coupled perturbative description 
in a gravity theory. The gauge/gravity duality provides the possibility for constructing a 
calculable theory of the strong interactions at low energies. Understanding the nonpertur- 
bative low-energy physics of quantum chromodynamics (QCD) is crucial, for instance, to 
formulating a quantitative description of the mass spectra of bound states of quarks and 
gluons. 

We proposed in [4, 5J a gravity dual to pure M = 1 SU {N) gauge theory in four dimensions 
in the large limit in the IR which reproduces the renormalization group fiow and the 
pattern of chiral symmetry breaking of the gauge theory and produces linear confinement of 
quarks within relativistic quantum theory of ten-dimensional superstrings for the first time. 

In this note we study fiuctuations of the dilaton on the gravity background to calculate 
mass spectrum of J^*" = 0+"*" glueballs, spin-0 bound states of gluons which are even under 
parity and charge conjugation. The masses are given by 

mn = - A„, n = l,2, (1.1) 

27rr^ 

where |^ is a combination of the parameters of the theory and are the roots of the 
function 

/(A) = M.^,,_.(A) + ^ M.,_.(A) , (1.2) 
where Ma , , _i (A) and M\ _i (A) are Whittaker functions of the first kind. 

4"""-^' 4 ^ ' 4 ' 4 ^ 
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The spectrum is discrete. The ratios of the masses are independent of the parameters 
and serve as predictions for 0^^ glueball masses in A/" = 1 SU (N) gauge theory in the large 
N hmit. Comparison of the ratios to recent large non-supersymmetric lattice QCD data 
available for the lowest three states shows agreement to within five percent. 



2 Supergravity background 

In this section we briefly summarize features of the supergravity background in [U |5]. The 
starting point of the construction is the supergravity background of unwrapped and fiat 
D7-branes on C^/Z2 with running axion and dilaton fields and the eight-dimensional 
A/" = 1 SU{N) gauge theory that lives on the branes. The background geometry of the 
eight-dimensional gauge theory is noncompact with a conic singularity. The D7-branes 
are then wrapped over a 4-cycle on <C^/Z2 x x T^/Z2 to construct A/" = 1 SU{N) 

gauge theory in four dimensions. The wrapped D7-branes turn on all Fi, F3, H3, and F5 
fluxes of type IIB theory and induce torsion. A/" = 1 supersymmetry in four dimensions is 
preserved and the background is stable by a balance between fluxes and torsion. The conic 
singularity in the background of the eight-dimensional gauge theory is smoothed out by the 
fluxes in the background of the four-dimensional gauge theory on the wrapped branes. The 
supergravity solutions are obtained by solving the supersymmetry equations in [B] which were 
generalizations of the equations in [H |S] . The solutions to the supersymmetry equations with 
the Bianchi identities imposed on them then automatically solve the bosonic supergravity 
equations of motion. 

The extra six-dimensional space of the ten-dimensional spacetime is parameterized by 
one radial coordinate r with range < r < r^e^ss^ and five angles ip, (pi, ip2, ^ps, and ip^, 
each having the range between and 27r; r and ip are coordinates on C^/Z2 ~ x and 
the remaining angles are coordinates on T'^/Z2 x T^/Z2. We call r = the IR boundary 
and r = rsC^asN ^j^e UV boundary; QgN is the 't Hooft coupling at the UV boundary. The 

D7-branes are wrapped over T^/Z2 x T^/Z2 at r = r^. The background metric is 

(is°io^ = cosh u dx\ 3 + sech u {dp^ + d^''^ 

+dipl + d^l + d^l + d^l)^ , (2.1) 

where 



sechM= Wl-(^p)2, p = 31n(-), (2.2) 
and dx\^ is the metric on fiat four-dimensional spacetime R^'^ with coordinate x^. We use 
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uppercase indices M, N, - ■ ■ to denote the coordinates of the ten-dimensional spacetime. A 
superscript of is put to denote background metric and fluxes excluding fluctuations. 
The dilaton and the background fluxes are 



where 



F3 = 1^ tanh u {dip A d(pi A d(fi2 + dip A dips A dip^) , 
Hl = 3^rl Vl + 2cosh2M {dt A dipi A dip2 + dt A dips A dcp^) , 
FO = (1 + ^10)^5° = i^r^ (1 + tanh^ v) dip A difi A dip2 A dtp^, A ^934 

Vl + 2 cosh 2m cosh 2m cosh^-ud^a; A dt, (2.3) 



dt^ , (2.4) 
Vl + 2 sech 2m 



The background metric and fluxes above solve all the bosonic supergravity equations 
given by 

dM (v^e-*G^^a^$) = -\fI - + l^e-^'Hl (2.5) 



Run = -4aM$a;v* + ^e2*aMCo9^Co + ^(i/3)MOp(i^3)jv'''' 

+ |e^*(-^3)MOp(-^3)Ar'^^ + -^e^^ {Fb) MOPQR{Fb) N^^^^ 

-G„.(^i^l + ^e»F| + 4e»F|), (2.6) 
(iFa = -Fi A t^^5 = ^3 A F3, (2.7) 



d ^10 (i^3 - e^^CoFa) = -e^^Fg A Fa. (2.8) 

The radial coordinate r is mapped to the scale A of the gauge theory as A/A^ ~ r/rs, 
where A^ is the nonperturbative scale at which the gauge coupling formally diverges. The 
coordinate ip is mapped to the Yang-Mills angle in the gauge theory as © = Nip. The 
gauge coupling of the of the four-dimensional gauge theory is related to the gauge coupling 
of the eight-dimensional gauge theory and the volume of the 4-cycle wrapped by the D7- 
branes, and it runs by inheriting the running of the dilaton from the background of the 
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eight-dimensional gauge theory. The axion potential and the dilaton are given by Co = 2:^1/^ 
and e~* = ^ p in the background of the eight-dimensional gauge theory. The background 
geometry of the four-dimensional gauge theory is compact and conformally Calabi-Yau. The 
't Hooft coupling of the four- dimensional gauge theory is given by = which is large 
in the IR and equals QgN = e*°A^ at the UV boundary. The curvature is smallest in the IR 
region where the gauge theory is strongly coupled and a dual gravity description is useful. 
The internal space normal to the wrapped D7-branes at the IR boundary is whose radius 
is set by the nonperturbative scale of the gauge theory independent of the 't Hooft coupling 
and spacetime is four- dimensional at the UV boundary. The potential energy of a quark- 
antiquark pair located at the UV boundary increases linearly with their spatial separation 
in four dimensions. 



3 0++ glueball masses 

We want to study 0"*"^ glueballs in four dimensions. The masses of 0"^"^ glueballs in the 
gauge theory are related to the two point correlation function of the square of the gluon field 
strength tensor (trF^trF^). Following a prescription of the gauge/gravity duality PQ 121 E], 
the glueball states are mapped to supergravity states that couple to the operator trF^. 
Because the dilaton couples to trF^, the 0+"^ glueball masses can be computed by studying 
linear fluctuations of the dilaton [21 HO]- To write appropriate equation of motion, we start 
with $ = $0 + 5$, Fi = FO + 6F1, F3 = + H3 = Hi + 6Hs, F5 = + and 
Gmn = G\.^^ + 6Gmn, where the second terms are the fluctuations. Because the dilaton 
field in the background of the four-dimensional gauge theory is constant, e*° = gs, the left 



hand side of (2.5) contains only fluctuations of the dilaton but not of the metric at linear 



order and we write (2.5) as 

dM{\^ G^'^'^Sjv^*) = 0, (3.1) 

where linear order fluctuations in the fluxes and the metric need to satisfy 1^2(12 e'^'^" 6{F^) + 

g2<i> (5(^1) — = (5$ in order for the background with the linear fluctuations to pre- 

serve supersymmetry. The equations for the fluctuations of the metric and the fluxes would 
then involve 5$. Our interest here is to study fluctuations of the dilaton governed by the 



decoupled equation given by (3.1) 



To find the mass measured by a four-dimensional observer, we make a plane wave expan- 
sion of the form 6^{x,p) = (f){p) e^^'^ and the glueball mass m is given by = —k'^. Using 
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the metric given by (2.1) in (3.1), 



<t>"{y) + X\l-y')m = ^. 



(3.2) 



where 



2'Krs 
9sN 



gsN 



and < y < 1. But (3.2) can be written as a Whittaker equation 



(3.3) 



w"(z)+ -- + — + 



4 4z IQz^ 



w{z) = 0, 



(3.4) 



where 



(3.5) 



The general solution to (3.2) can then be written as 



\ 4 ' 4 4 ' 4 



(3.6) 



where Ma i (Ay^) and + Wx i (Xy"^) are respectively Whittaker functions of the first kind and 

414 414 

the second kind and C and D are constants. 

Now we need to impose appropriate boundary conditions at y = and y = 1. We demand 
that physically acceptable modes of fluctuation be smooth at both boundaries, which requires 
Neumann boundary condition at both ends, (f)'{0) = 0'(1) = 0. 

First, the boundary conditions accommodate a massless mode with constant profile, 
which can also be read off directly from (3.2), (j)o{y) = c, where c is a constant. Because the 
classical U{1) R-symmetry of the gauge theory is anomalous and reduced to a discrete Z2N 
symmetry in the quantum theory, its breaking to Z2 by gaugino condensation does not lead 
to a massless Goldstone boson. Therefore, we avoid this mode by setting c = 0, since it is 
constant and nondynamical. To find the massive modes, first we impose 0'(1) = and (3.6) 
reduces to 



c 



(A + 3)Ma^i i(A) + (A - 1)Ma i(A) 



AWx 

4 



iX) - {X - l)Wx ,{X) 

4 ' 4 



Wx 1 

4 : 4 



iXy' 



(3.7) 



We then impose 0'(O) = to obtain discrete values of A„ as roots of the function /(A) 
riven by (1.2), and the corresponding modes (j)n{y) are obtained using A„ for A in ( 3.7[ ). The 
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(a) n = 1. 




(c) n = 3. 




(b) n = 2. 




(d) n = 4. 



Figure 1: The first few modes (j)n{y) plotted using the expression given by (3.7) with the 



roots A„ obtained from (1.2) for A. 



first few numerical values of A„ are given in Table [T] The roots are approximately given by 



Ai 


A2 


A3 


A4 


As 


As 


A7 


As 


A9 


Aio 


4.287 


8.304 


12.311 


16.315 


20.317 


24.319 


28.320 


32.321 


36.322 


40.323 



Table 1: The first few roots of /(A). The glueball masses are given by (1.1). 



4.3 + 4(n -l),n = 1,2, 



One way to obtain (1.2) is to make a series expansion 



of (3.7) in y, set the coefficient of the linear term to zero, and use the recursion relation 
W^fc,/.(A) = r(T/^zj!z^Mk,^{\) + r(i/2+^Lfc) ^fc,-M(A) that is vahd for non-integer 2fi, which is 
the case here. The first few modes, normalized as dy (j)n{yY = plotted in Figure [l] 



using (3.7) with the numerical values for A„ obtained from the roots of /(A). 

The glueball wavefunctions are nonzero and have largest magnitudes, and increasingly 
so for larger n, at the location of a four-dimensional observer at the UV boundary. On the 
other hand, modes of fiuctuations in the angular directions near the UV boundary, where the 
radius of the transverse internal space gets smaller and spacetime becomes four- dimensional, 
are much heavier. 
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Figure 2: The first few Vn{y)- Tlie glueball masses are given by (1.1), wfiere = V^(l) are 



tlie roots of (1.2). 



The equation of motion of (f){y) given by (3.2 ) can be interpreted as a Shrodinger equation 
of the form (f)"{y) - XV{y)(f){y) + XE(j){y) = with V{y) = Xy^ and E = X = V{1). The 
different values of A„ can then be viewed as energies corresponding to different modes (f)n{y) 
for potentials Vn{y) with depths which increase with increasing A„. Plots of the first few 
Vn{y) are shown in Figure |2] 



4 Comparison to lattice QCD data 



The ratios of the masses are independent of the parameters of the theory, = ^r^, 

and serve as predictions. A lattice computation for the lowest three O"*"^ glueball masses in 
non-supersymmetric SU{oo) QCD is given in [H], with the states denoted by A^^, A^^*, 
and A'f'^**. This allows us to compare two independent mass ratios as shown in Table [2j 
Both ratios agree to within five percent. 





Supergravity 


Lattice 


mi 


8.304 _ qQ7 
4.287 -L--''-" 


1-456 _ 1 071 
0.778 -L-O'-L 


ms 

7712 


M = 1-483 


1.456 -L-'iJ-u 



Table 2: Comparison of mass ratios between our supergravity theory and the non- 
supersymmetric SU {00) lattice QCD data in [llj . 
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5 Conclusions 



The supergravity background that we argued in [4' ;5J to correspond to pure N = 1 SU{N) 
gauge theory in four dimensions in the large hmit in the IR produces mass spectrum of 
0"'""'" glueballs that is quahtatively consistent with expectations from the gauge theory side 
and quantitatively in striking agreement with available data from non-supersymmetric large 
N lattice QCD. 

Because the background geometry is compact, the discreteness of the spectrum follows 
with appropriate boundary conditions imposed at the IR and the UV boundaries. The 
magnitudes of the masses become smaller as the 't Hooft coupling at the UV boundary gets 
smaller, consistent with expectation from the gauge theory side where a smaller 't Hooft 
coupling corresponds to smaller effective gauge coupling in the large limit. From the 
gravity theory point of view, a smaller 't Hooft coupling at the UV boundary corresponds to 
a larger size of the compact background and smaller curvature which gives modes with smaller 
masses. A four dimensional observer at the UV boundary measures a 't Hooft coupling of 
magnitude QsN that can be chosen appropriately. Furthermore, Kaluza-Klein excitations 
in the angular directions near the UV boundary are much heavier and a four-dimensional 
observer sees only the radial excitations. 

The ratios of the masses are independent of the parameters of the theory and the interval 
between consecutive masses is nearly constant and approximately 4(|^) which can be cho- 
sen appropriately to produce absolute magnitudes of the masses. Comparison of the mass 
ratios to recent non-supersymmetric large A^ lattice QCD data available in |llj for the first 
three masses shows agreement of the two independent ratios to within five percent. The su- 
pergravity expressions for the masses and wavefunctions in this note include the higher mode 
0"'""'" glueballs and additional lattice data would be quite interesting for further comparison. 

It would also be interesting to study metric and flux fluctuations on the background and 
calculate the spectrum of other J^*^ glueballs and to investigate their relations to the O^"*" 
glueballs and see if they fall on Regge trajectories. 

We would like to re-emphasize that the supergravity background presented in ^ |5j and 
briefly summarized in Section |2] has a number of remarkable features. It reproduces the 
renormalization group flow and the pattern of chiral symmetry breaking of the gauge theory. 
It produces linear confinement of quarks in four dimensions within relativistic quantum 
theory of ten-dimensional superstrings. The background is compact and the curvature is 
smallest in the IR where the gauge theory is strongly coupled and a dual gravity description 
is useful. The radius of the transverse internal space at the IR boundary is set by the 
nonperturbative scale of the gauge theory and is independent of the 't Hooft coupling. 
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Smaller 't Hooft coupling at the UV boundary corresponds to smaller curvature of the 
background. The background preserves supersymmetry and is stable. 

The linear confinement of quarks the background produces together with the agreement of 
the ratios of the 0'^''" glueball masses obtained in this note with available non-supersymmetric 
large N lattice QCD data means to us that J\f — 1 SU (N) gauge theory has similar mass 
spectra as in QCD and the supergravity theory can be used for exploring additional nonper- 
turbative phenomena in QCD. 
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